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It has recently been shown that a Hagedorn phase of string gas cosmology can provide a causal 
mechanism for generating a nearly scale-invariant spectrum of scalar metric fluctuations, without 
the need for an intervening period of de Sitter expansion. In this paper we compute the spectrum 
of tensor metric fluctuations (gravitational waves) in this scenario, and show that it is also nearly 
scale-invariant. However, whereas the spectrum of scalar modes has a small red-tilt, the spectrum 
of tensor modes has a small blue tilt, unlike what occurs in slow-roll inflation. This provides a 
possible observational way to distinguish between our cosmological scenario and conventional slow- 
roll inflation. 
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Introduction. Strinff gas cosmology (SGC) 0,|3| (see 
also I3 for early work, |j,|5j for recent overviews, and [g 
for a critical review) is an early approach to string cos- 
mology, based on adding minimal but crucial inputs from 
string theory, namely new degrees of freedom - string os- 
cillatory and winding modes - and new symmetries - T- 
duality - to the now standard hypothesis of a hot and 
small early universe. One key aspect of SGC is that the 
temperature cannot exceed a limiting temperature, the 
Hagedorn temperature Th |3- This immediately pro- 
vides a qualitative reason which leads us to expect that 
string theory can resolve cosmological singularities [l|. 
From the equations of motion of string gas cosmology 
y, I3 it in fact follows that if we follow the radiation- 
dominated Friedmann-Robertson- Walker (FRW) phase 
of standard cosmology into the past, a smooth transi- 
tion to a quasi-static Hagedorn phase will occur. In this 
phase, the string metric is quasi-static while the dila- 
ton is time-dependent. Reversing the time direction in 
this argument, we can set up the following new cosmo- 
logical scenario [jj]: the universe starts out in a quasi- 
static Hagedorn phase during which thermal equilibrium 
can be established over a large scale (a scale sufficiently 
large for our current universe to grow out of it following 
the usual non-inflationary cosmological dynamics). The 
quasi-static phase, however, is not a stable fixed point of 
the dynamics, and eventually a smooth transition to a 
radiation-dominated FRW phase will occur, after which 
point the universe evolves as in standard cosmology. 

Recently |9|, llO| it was discovered that string thermo- 
dynamic fluctuations during the Hagedorn phase lead to 
scalar metric fluctuations which are adiabatic and nearly 
scale-invariant at late times, thus providing a simple al- 
ternative to slow-roll inflation for establishing such per- 
turbations. It is to be emphasized that this mechanism 
for the generation of the primordial perturbations is in- 
trinsically stringy - particle thermodynamic fluctuations 
would lead to a spectrum with a large and phenomeno- 
logically unacceptable blue tilt. 



We briefly remind the reader of the key features of 
the new structure formation scenario. The Figure shows 
a sketch of the space-time evolution: at early times t 
(t < tji) the universe is in the quasi-static Hagedorn 
phase. The physical wavelength of any perturbation 
mode (characterized by having constant momentum k in 
comoving coordinates) is approximately constant. The 
key point is that the Hubble radius (which sets the limit 
on which causal processes can locally set up fluctuations 
- see e.g. [ll| for a concise overview of the theory of cos- 
mological perturbations and [l3| for a comprehensive re- 
view) is essentially infinite, thus allowing a causal mecha- 
nism for the generation of primordial fluctuations. Near 
t = tji, a smooth transition from the Hagedorn phase 
to the radiation-dominated phase of standard cosmology 
occurs. The Hubble radius decreases rapidly to take on 
a minimal value which is microscopic (set by the temper- 
ature at t = tji which will be close to the Hagedorn tem- 
perature). Thus, fluctuation modes of relevance to cur- 
rent cosmological observations exit the Hubble radius at 
times ti{k) close to tji. During the radiation-dominated 
FRW phase, the Hubble radius increases linearly in t, 
while the physical wavelength of a perturbation mode 
grows only as i^^^. Thus, at late times tf{k) the fluctua- 
tion modes will re-enter the Hubble radius. Since the pri- 
mordial perturbations in our scenario are of thermal ori- 
gin (and there are no non- vanishing chemical potentials), 
they will be adiabatic, and since they propagate on super- 
Hubble scales for a long time during the FRW phase, they 
will be squeezed and will lead to the same type of acous- 
tic oscillations in the angular power spectrum of the cos- 
mic microwave background (CMB) anisotropics as what 
is produced in slow-roll inflation models. 

In this Letter, we generalize our previous analysis 
[si UJI to allow us to compute, in addition to the 
scalar metric fluctuations (the metric perturbation modes 
which transform as scalars under spatial rotations and 
which couple to matter), also the tensor modes (gravi- 
tational waves). We flnd that the resulting spectrum of 



tensor modes is also nearly scale- invariant, but that it 
has a slightly blue tilt, unlike what happens in slow-roll 
inflation where the tilt for the gravitational wave spec- 
trum is also red. The scalar to tensor ratio is calculable 
from the detailed dynamics of the system (we postpone 
this calculation to a foUowup paper). It is set by how 
close the temperature is to the Hagedorn temperature 
when the scales which are measured today exit the Hub- 
ble radius at the end of the quasi-static Hagedorn phase. 
It is not suppressed by the same small parameter which 
sets the amplitude of the spectrum of scalar metric fluc- 
tuations. Thus, a measurement of the background of 
stochastic gravitational waves would allow us to distin- 
guish our scenario from the usual slow-roll inflationary 
models. 

The outline of this Letter is as follows: we begin with 
a quick recap of the method of [sl lHI to compute mass 
fluctuations in the Hagedorn phase of SGC and general- 
ize the method to yield fluctuations of arbitrary compo- 
nents of the energy-momentum tensor of the string gas. 
Then, we relate the spectrum of gravitational waves in 
late time cosmology to the fluctuations of the energy- 
momentum tensor in the Hagedorn phase and show that 
the resulting power spectrum is nearly scale-invariant, 
like the spectrum of the scalar metric fluctuations. 

In the following, we assume that our three spatial di- 
mensions are already large during the Hagedorn phase 
(for a possible mechanism to achieve this see PJi), while 
the extra spatial dimensions are confined to the string 
scale. For a mechanism to achieve this in the context of 
SGC see [l| (see however Il4l,'l5i| for caveats), and for a 
natural dynamical mechanism arising from SGC to sta- 
bilize all of the moduli associated with the extra spatial 
dimensions see [3, E3, El El 113 ■ To be specific, we 
take our three dimensions to be toroidal. The existence 
of one cycles results in the stability of string winding 
modes - and this is a key ingredient in our calculations. 

Energy-Momentum Tensor Correlation Func- 
tions for Closed Strings. We begin with the thermal 
canonical partition function Z of a gas of closed strings, 
from which the mean energy-momentum tensor {T^^) 
can be found by 



{T^.) = 2 



G^" d\nZ 



Vi'lnZ . 



(1) 



where Gfj,i, is the Euclidean metric of space-time (the 
time coordinate is compactified to a circle of radius /3, 
where /3 is the inverse temperature. 

Our aim is to calculate the fiuctuations of the energy- 
momentum tensor on various length scales R. For each 
value of R, we will take our spatial coordinates to run 
over a fixed interval, e.g. [0, 27r]. Thus, the metric will 
have dimensions and will be given by 



G^, = diag[/32, i?2^i?2,i?2] 



(2) 



Letting another derivative operator V act on In Z gives 
two terms, one for which both derivative operators act 
on Z, a second which will contain the product of terms 
where one V acts on Z. Specifically (and symmetrizing 
over the indices), we find the mean square fluctuation to 
be 

G^% = (<5r^,5r^) = (r^r^)-(r^)(r^) 




dlnZ 



G^* d\nZ 



GdG^'^ 



(3) 



with(5T^, =r^, -(T^,). 

The partition function Z = exp {—(3F) is given by the 
string free energy F. Thus, the string thermodynamical 
fluctuation in the energy density, here denoted by the 
correlation function G^o^'oi becomes 



C%% = {Sp') = {p') {pf 



R^dp 



F + (3 



dF 

'dp 



i?6 



G 



V ■ 



(4) 



where Gy = {d{E)/dT) with E = F + P{dF/d(3) and 
V = R^ is the volume of three compact but large spatial 
dimensions. The fluctuation in spatial diagonal parts of 
the energy-momentum tensor can be found by 



i i 



(ST' 

1 



(T* 



d 



1 



(T^)2 
dF 



1 dp 



(5) 



/?i?3 a In i? V R^d InRj I3R'^ OR 
with no summation on i. The pressure p is given by 

p = -~V-^{dF/d\nR) = T{dS/dV)E. (6) 

In the following, we will compute the two correlation 
functions Q and jSJ) using tools from string statistical 
mechanics. Specifically, we will be following the dis- 
cussion in J23 (see also [H H IH ill). The starting 
point is the formula S{E,R) — lnn{E,R) for the en- 
tropy in terms of n{E,R), the density of states. The 
density of states of a gas of closed strings on a large 
three-dimensional torus (with the radii of all internal di- 
mensions at the string scale) was calculated in 22] and 
is given by 

n{E, R) ~ /3ff e'3-^+"-^ [1 + 517(1) {E, R)] , (7) 

where 551(i) comes from the contribution to the density 
of states (when writing the density of states as a Laplace 
transform of Z{(3), which involves integration over /3) 
from the closest singularity point (ii to (in = {^/Th) 
in the complex fi plane. Note that /3i < Ph, and /3i is 
real. From |23| we have 



5^[i)[E,R) 



5! 



(8) 



In the above, n^j is a (constant) number density of order 
i~^ {£s being the string length) and pH is the 'Hagedorn 
Energy density' of the order i^^, and 



(3h-Pi 



{i^/R'^), fori?>£,, 
(Ryes), fori?«4. 



(9) 



To ensure the vahdity of Eq. ((Tj) we demand that S^(^i^ <C 
1 by assuming p = (E/V) » pH- 

Combining the above results, we find that he entropy 
of the string gas in the Hagedorn phase is given by 

S{E, R) ~ (3hE + nnV + In [l + (50(1,] , (10) 

and therefore the temperature T{E,R) = [{dS/dE)vY'^ 
will be 



T: 






T, 



H 



Ph-P- 



Ph 



-5^,1 



(1) 



(11) 
Using this relation, we can express 5D,ii\ in terms of T 
and R via 



tl6n.,^-^(i-^ 



Th 



In addition, we find 

{E) ~ l^^R^ In 



Th 



(.It 



R^I-T/Th) 



(12) 



(13) 



Note that to ensure that 551(1) <C 1 and (E) :^ PhR^, 
one should demand (1 - T/Th) > (i'i/R^)- 

The results IJIU]) and ifT^ now allow us to compute the 
correlation functions Q and |SJ- We first compute the 
energy correlation function Q), noting that 



Cv 



rpZ 



d^S{E,R) 
dE^ 



i?V^3 



T (1 - T/Th) ' 

(14) 



from which we get 

C%% = {Sp') 



£lil-T/TH)R'' 



(15) 



Note that the factor (1 — T/Th) in the denominator is 
responsible for giving the spectrum a slight red tilt. It 
comes from the differentiation with respect to T. 
Evaluating © 



fj^ m^ T^ 2 (1 - T/Th) , 
p{E, R) w uhTh - j^ In 



^T 



R^1-T/Th) 



immediately yields 

r(i - t/Th) 



c 



^3i?4 



In^ 



(1 - T/Th) 



(16) 



(17) 



Note that since no temperature derivative is taken, the 
factor (1 — T/Th) remains in the numerator, and it is this 



fact which will lead to the slight blue tilt of the spectrum 
of gravitational waves. 

Tensor Modes from Hagedorn Fluctuations. In 

this section we estimate the dimensionless power spec- 
trum of gravitational waves. First, we make some gen- 
eral comments. In slow-roll inflation, to leading order 
in perturbation theory matter fluctuations do not couple 
of tensor modes. This is due to the fact that the mat- 
ter background field is slowly evolving in time and the 
leading order gravitational fluctuations are linear in the 
matter fluctuations. In our case, the background is not 
evolving (at least at the level of our computations), and 
hence the dominant metric fluctuations are quadratic in 
the matter field fluctuations. At this level, matter fluctu- 
ations induce both scalar and tensor metric fluctuations. 
Based on this consideration we expect that in our string 
gas cosmology scenario, the ratio of tensor to scalar met- 
ric fluctuations will be larger than in simple slow-roll in- 
flationary models. 

We will extract the amplitude of the gravitational 
wave spectrum from the spatial fluctuations C*j*j of the 
energy-momentum tensor. Strictly speaking, it is the off- 
diagonal components which will couple uniquely to the 
tensor modes. We will estimate their order of magni- 
tude by the order of magnitude of the diagonal terms 
computed in the previous section. This gives a good ap- 
proximation, as can be checked by letting the metric in 
^ depend on three separate scales Ri (where the index i 
runs from 1 to 3), and by extracting the off-diagonal cor- 
relation functions following the method of the previous 
section, but taking mixed spatial derivatives. 

Tensor perturbations in a spatially flat FRW universe 
take the form 



ds^ = -dt^ + a^{t)(5ij + hij)dx^dx^ 



(18) 



Since to linear order in hij the Einstein tensor for fluctu- 
ations on a scale k is proportional to k'^hij{k), it follows 
from the space-space Einstein equations that 



K ii"ij \'^) 



(87rG)(5Ty(fc). 



(19) 



The power spectrum of the right-hand side of the above 
equation is given by the correlation function C^/j- Thus, 
the dimensionless gravitational wave power spectrum is 
given by this correlation function. Therefore, from H19() 
one can calculate the dimensionless power spectrum for 
h^ , where h^ is the amplitude of either of the two grav- 
itational wave polarization modes. Dropping the super- 
script ± (due to symmetry both polarization modes will 
be equally excited) we obtain 

fc3|;i(fc)|2 ^ k-\8TrG)^C'j'j. (20) 

Inserting our result ^ for the correlation function yields 



P|/i(/c)P 



GAtt'^G'^T 



(I-T/Th) \n 



£W 



(1 - T/Th) 



(21) 



which, for temperatures close to the Hagedorn value re- 
duces to 



fc3|/i(/c)|^ 



'^pi 



{l-T/TH)\n 



1 



PM- 



■{1-T/Th) 



(22) 

This shows that the spectrum of tensor modes is - to a 
first approximation - scale-invariant. 

Discussion. Our result (|22|l for the power spectrum 
of gravitational waves should be compared to the result 
for the power spectrum of scalar metric fluctuations com- 
puted in [21 : 



V^{k) 



L 



1 



(23) 



1 - T/Th ■ 

Note that for a fixed scale fc, both ^2'A\ and 1)22(1 must be 
evaluated at the time ti{k) when the mode k exits the 
Hubble radius at the end of the Hagedorn phase. Since 
ti{k) increases slightly as k increases, the temperature 
T{ti{k)) will be slowly decreasing. Hence, the expression 
in front of the logarithm in our final expression H22|) for 
the power spectrum of tensor fluctuations yields a slight 
blue tilt. Values of k for which the perturbative analysis 
of string gas cosmology is consistent are on the high k 
side of the zero of the logarithm in 1)22(1 - this follows from 
the condition |(551(i)| ^ 1 and ((12(1 . Hence, the logarithm 
factor in 1(22(1 adds to the blue tilt of the spectrum. 

A heuristic way of understanding the origin of the 
slight blue tilt in the spectrum of tensor modes is as fol- 
lows. The closer we get to the Hagedorn temperature, 
the more the thermal bath is dominated by long string 
states, and thus the smaller the pressure will be com- 
pared to the pressure of a pure radiation bath. Since the 
pressure terms (strictly speaking the anisotropic pressure 
terms) in the energy-momentum tensor are responsible 
for the tensor modes, we conclude that the smaller the 
value of the wavenumber k (and thus the higher the tem- 
perature T{ti{k)) when the mode exits the Hubble radius, 
the lower the amplitude of the tensor modes. In contrast, 
the scalar modes are determined by the energy density, 
which increases at Titiik)) as k decreases, leading to a 
slight red tilt. 

It is also interesting to consider the tensor to scalar 
ratio r. Comparing ((22(1 and ((23(1 we see that this ratio, 
evaluated on a scale k is given by 



r ~ {l-T{U{k))lTHf\n 



1 



£2fc2 



(1 - T{U{k))/TH) 



In principle (if the dynamical evolution from the Hage- 
dorn phase to the radiation-dominated FRW phase were 
under complete analytical control) this quantity is cal- 
culable. If the string length were known, the factor 
(1 — T/Th) could be determined from the normaliza- 
tion of the power spectrum of scalar metric fluctuations. 
Since the string length is expected to be a couple of or- 
ders larger than the Planck length, the above factor does 



not need to be very small. Thus, generically we seem to 
predict a ratio r larger than in simple roll inflationary 
models. 

Based on the results of this Letter it thus appears 
promising that our scenario will give rise to testable pre- 
dictions. 
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